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Abstract 

Linked cluster expansions provide a useful tool for both analytical and 
numerical investigations of lattice field theories. The expansion parameter(s) 
being the interaction strength(s) fields at neighboured lattice sites are coupled, 
they result into convergent hopping parameter like series for free energies, cor- 
relation functions and in particular susceptibilities. We consider scalar fields 
with 0(N) symmetric nearest neighbour interactions on hypercubic lattices 
with possibly finite extension in some directions, thus including field theories 
at finite temperature T. We improve known and develop new techniques and 
algorithms to increase the order n the expansions can be computed to in such a 
way that detailed information on critical behaviour can be extracted from the 
susceptibility series. This concerns both simple moments as well as higher cor- 
relations such as 4- and 6-point functions used to define renormalized coupling 
constants. Particular emphasis is done on finite temperature field theory. In 
order to be able to measure finite temperature critical behaviour, the order of 
explicit computation n has to be sufficiently large compared to in lattice 
units. 2- and 4-point susceptibility series are computed up to and including 
the 18th order and beyond. 
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1 Introduction 



Most physical problems do not allow for closed analytic solutions. This is mainly 
due to the fact that realistic models can and do evolve considerable complexity. 
This is particularly the case with statistical mechanics and field theory, where the 
number of degrees of freedom becomes large or even infinity. Close to criticality 
those degrees of freedom interact collectively. One is forced to use approximate 
methods. 

A desirable practice is to construct and apply systematic expansions about ex- 
actly solvable cases such as harmonic oszillators or zero-dimensional field models. 
There are essentially two categories. On the one hand, we use asymptotic methods. 
Common examples are saddle point expansions and weak coupling perturbation 
theory. Those methods are designated to an accurate quantitative description of 
particular sectors, such as the short distance behaviour of QCD. Their fait is that 
they are only asymptotically convergent. On the other hand, if available, conver- 
gent expansions have the advantage that in principle they allow for arbitrary precise 
measurements within the domain of convergence. The price to be paid is that those 
expansions typically are more involved, need their own techniques, and accurate 
results actually need high orders to be carried out. 

We are concerned here with the technique of linked cluster expansion (LCE) 
applied to scalar fields on a hypercubic lattice. It is a special form of the hopping 
parameter expansion for the computation of the free energy, connected correlation 
functions and susceptibilities. The expansion is a power series in the strength(s) 
fields at different lattice sites are coupled. For standard actions this just is the 
nearest neighbour coupling. The high order behaviour of correlation functions and 
susceptibilities allows for a determination of critical couplings and critical exponents. 

A good deal to manage the complexity of extreme orders is to use graph theo- 
retical methods. To every finite order a correlation function is represented as a sum 
over a finite number of graphs with appropriate weights depending on the geometry 
of the lattice and the couplings involved. Calculation of extreme orders need tech- 
niques and formalization so that the explicit computation in the end can be done 
by a computer. So the main task is the developement of efficient algorithms and 
implementations. 

The techniques of LCE are of course not new. An elementary introduction can 
be found in [|l], 0. The algorithmic formulation of this expansion for general 0(N) 
symmetric scalar field models has been pioneered by Liischer and Weisz in the famous 
work 0, in the course of the solution of lattice $^-theory in four dimensions 0, |^. 
They succeeded to compute two- and four-point susceptibility series to 14th order 
in the hopping parameter on infinite lattices. 

Our main emphasize here is to further develop and improve the techniques and 
methods. The aim is to increase the order of explicit computation for field models 
both with discrete and continuous symmetry group in such a way that detailed 
knowledge on critical behaviour can be extracted from the LCE series. A major 
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field of application we have in mind are scalar field theories at finite temperature. 
Those models live on a lattice with toroidal symmetry in temperature direction 
with periodic boundary conditions imposed, of extension Lq given by the inverse 
temperature in lattice units, Lq = T~^. Graphs are able to "feel" temperature only 
if they can wind around the torus. There is a least order nmin depending on Lq, the 
onset of T-dcpcndcncc. We expect that njni^ has to be sufficiently large compared 
to Lq. That is, critical behaviour at finite T is postponed to higher orders compared 
to the case of a lattice infinitely extended in all direction. 

The organization of this article is as folUows. We start in Sect. 2 with a re- 
minder to the basics of the LCE of lattice scalar fields. This concerns motivation, 
ideas and the problems to be resolved. LCE series are represented as a sum over 
(equivalence classes of) appropriate graphs with associated weights. Definitions and 
graph theoretical fundamentals that are used throughout are introduced. Sect. 3 is 
dedicated to the algebraic representation of graphs by incidence matrices. This is 
a necessary business to be carried out in order to generate and compute diagrams 
numerically. The main problems to be fixed are uniqueness of the representation, 
i.e. 1-1 correspondence between equivalence classes of graphs and appropriately de- 
fined canonical incidence matrices, and the construction of very efficient algorithms 
to transform any incidence matrix into its canonical form. This is at the heart of 
the problem and necessitates to work around huge factorials. Introducing extended 
order relations on vertices is the way out. Sect. 4 describes the construction of all the 
graphs of interest, such as connected and IPI ones, along an inheritance tree from 
so-called "base classes" of less complexity and least number of graphs. The order 
by order generation is done within the simplest base class only. The computation 
of the weights belonging to a graph is discussed. Sect. 5 summarizes the results. 

2 Linked Cluster Expansion 
2.1 Notations and principles 

We consider a hypercubic D-dimensional lattice A, either equal to Aq = or to 
Ai = Z/Lo + Z^~\ where D, Lq G N and Lq is an even number (finite temperature 
lattice with temperature T = Lq^ in lattice units). Periodic boundary conditions are 
imposed. The restriction to even Lq is for convenience and leads to a considerable 
reduction of the number of graphs, as will be seen below. To be specific, the models 
we discuss are described by the partition function 

Z{J, v) = I H d^'<!?{x) exp (-5($, v) + Yl ■ *(^))' (1) 
where $ denote a real, A'"-component scalar field, J are external sources, and 

AT 
a=l 
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The action is given by 

1 ^ 

X ^ XytyeAa,b=l 

o 

The "ultra- local" (0-dimensional) action S should be 0{N) invariant and respect 
stability of the partition function (0) for sufficiently small Vab- For a $^ + $^ inter- 
acting theory, the standard canonical form is given by 

^($) = $2^Ai(<l'2-l)2 + A2($'-l)', (3) 

with A2 > or A2 = 0, Ai > 0. Dependence on the bare ultra-local coupling constants 
A is suppressed in the following. Fields at different lattice sites are coupled by the 
hopping term v{x,y). For the case of nearest neighbour interactions, 

, s \ 2k 5ab , X, y nearest neighbour, , 
"'^^^^'^^ = 1 '0 , otherwise. 

5 is the Kronecker symbol. On Ai, the nearest neighbour property is modulo the 
torus length. 

The generating functional of connected correlation functions (free energy) is given 

by 

W{J,v) = \nZ{J,v), (5) 

W'if^.L (a^l. • • • . ^2n) = < (Xi) ■ ■ ■ (X2„) 

and by a Legendre transform 

T{M,v) = W{J,v) -J2 Ji^) ■ M{x) 

= E rifi^Jxi,...,X2„)M,,(xO---M,,„(x2„), (6) 

n>0 ai,...,a2„ 

Ma{x) = , a = l,...,n 

dJa{x) 

the generating functional of 1-particle irreducible (IPI) correlations (vertex func- 
tional) is defined. In field theory they are used to define physical coupling constants, 
e.g. by 

ri?(p,-p) = -^{ml+p' + 0{p''))6a,b asp = (po = 0,p-.0), (7) 
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where denotes Fourier transform, and 



ritL(0, 0, 0, 0) = -J^^C4(a,6,c,d), (8) 

— —( 



rSU(0,---,0) = -^^C,{a^,...,a,). (9) 



The C„ denote the 0(N)-invariant tensors, 

C2(a,b) = 5a,b, (10) 

2n 

C2n('2l) • • • ) "^271) = X/ '-'2n-2(fl25 • • • 5 Oj, . . . , a2n), "H- > 2, 



i=2 



(^) denotes omission. The more directly accessible objects from the point of view 
of the linked cluster expansion are connected correlation and their susceptibilities. 
The latter are defined by 

Sa,bX2 = tyi?(0,0)=^<$„(x)$,(0)>^ (11) 



5a,, 1^2 = E^(^) < M^)MO) >'^= -Dr W^^p, -p)\^^^ , (12) 



dpi 

with g{x) = Ej=o Dr ^ D on Aq and g{x) = J2f=i xj, Dr ^ D - 1 on Ai, 

^C4(a,6,c,rf)x4 = W^a6cd(4)(0, 0,0,0) (13) 
^Ce(a,,...,ae)xe = W^i?..„e(0, ■ ■ ■ , 0). (14) 
The physical coupling constants as defined above are related to them by 

(15) 
(16) 

X4, (17) 

(xe - log . (18) 

If the non-ultralocal part of the action is switched off, i.e. v — 0, S{^,v = 0) = 
J2x S {^{x)), in turn W(J,v ^ 0) = Ex {J{x))- In particular, 

li^(2n) / ™ N _ J (2n-DM <^2n(ai, ■ ■ ■ , a2n) , for Xi = X2 = • • • = X2„, 

ai...a2„\-^i^T ■ ■ ■'•^2nJ ^ \ /■■ ,1 . 

, otherwise, 

(19) 





2Dr^, 
A*2 




2L>^^, 

/^2 




f2Dr 
V 


= 


f2Dr 
V /^2 
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with 

oc 



Q2n 



^2n= ITr^ W (J) 



(20) 

J=0 



The hnked cluster expansion is the Taylor expansion with respect to v{x,y) about 
this decoupled case, 



(21) 

V=0 

and the obvious generalizations to connected correlation functions according to (^. 
Multiple derivatives of W with respect to v{x,y) are computed by the generating 
equation 

d 1 / d^W dW dW \ 



dvab{x,y) 2 \dJaix)dJbiy) dJ^ix) dJb{y) 

The expansion ([2T| ) is convergent under the condition that the interaction v{x, y) 
is sufficiently local and weak. In finite volume this is easily shown, e.g. by means of 
the dominated convergence theorem. In order to get information on collective be- 
haviour, convergence has to be shown in the infinite volume case. To my knowledge, 
for this particular expansion the proof of uniform convergence for large volume has 
only recently been done P]. In many cases under consideration one restricts at- 
tention to couplings of the form (^, i.e to euclidean symmetric nearest neighbour 
interaction of strength k, only. By inspection, for the susceptibilities (|llD-(0) the 



coefficients of the series with respect to k are of equal sign. This remarkable prop- 
erty implies that the singularity closest to the origin lies on the positive real K-axis. 
This is identified with the critical point. In turn it is possible to extract quantitative 
information on behaviour even close to criticality from the high order coefficients of 
the series 0. 

The evaluation of the Taylor expansion rapidly becomes complex and cumber- 
some. A graph theoretical device is the appropriate tool to keep control of high 
orders. Correlation functions are represented as a sum over graphs with appropriate 
weights. We introduce some basic concepts. 



2.2 Some Graph Theory 

The notion of a graph used here is a slight specialization of the general definition 
[|, H adjusted for the LCE of susceptibilities of scalar field corrrelations. A graph 
or diagram is a structure 

T = iCr,Br,Er,^r), (23) 

where £r and Br ^ are disjoint sets, the internal lines and vertices of F, respec- 
tively. E-p is a map 

Er-.Br ^ {0,1,2,...}, 

V ^ Er{v), (24) 



6 



that assigns to every vertex the number of external hnes Et{v) attached to it. The 
number of external lines of F is given by J2veBr Er{v)- Finally, $r is the incidence 
relation that assigns internal hnes to their endpoints. More precisely, let P be the 
exchange of two vertices, 

P:BrxBr ^ Br x Br, 

{v,w) {w,v). (25) 

The imbedding map of ordered pairs of vertices onto unordered pairs is denoted by 
( " ), Br X Br = {Br x Br) /P. We identify {{v,v)\v e Br] with Br in the obvious 
way. Then 

^r-Cr^ {Br x Br) \ Br- (26) 

This maps every internal line onto its (unordered) pairs of endpoints. Lines are 
treated as unoriented. We notice the absence of tadpole lines, i.e. no hne has both 
of its endpoints equal. We furthermore set for every I e Cr with $r(0 — ('^^ '^) 

^r{l)={vM- (27) 

For completeness and later use we introduce further notions. External vertices 
are those that have external lines attached, 

Br,ext = {v&Br\Er{v)^Q}, (28) 

whereas internal vertices don't, i3r,mt = Br\ Br^ext- For every v,w E Br, 'm{v,w) 
denotes the number of common lines of v and w, i.e. the number of elements of 

M{v, w)^{le/:r\ ^r(0 = J^}- (29) 

The absense of tadpole lines imply that m{v, v) = for every vertex v. w is a 
neighbour of v if m{v, w) ^ 0, and 

Af{v) ^{w eBr\ m{v, w) ^ 0} (30) 

is the set of neighbours of v. It has N{v) elements. Multiple lines do not count here. 
For every integer n, n-vertices are those that have precisely n neighbours, 

B^'^ = {veBr\ N{v) = n}. (31) 

An n-vertex has at least n internal lines attached. 

Some topological notions and global properties will be of major interest in the 
following. A path is an ordered non-empty sequence p = {wi, . . . ,Wn) of vertices 
with m{wi, Wj+i) 7^ for alH = 1, . . . , n — 1. p is called a path in F from wi to Wn 
of length n — 1. If in addition Wn — wi, p is called a loop, with n — 1 its number of 
lines. If for every pair of vertices v,w & Br, v ^ w, there is a path from v to w, the 
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graph r is called connected. If every vertex v G Br has an even umber of internal 
and external lines attached, i.e. 



^('^) X! ''T^{'v,w) + Er{v) even, (32) 

the graph T is called even. In particular, then, F has an even number of external 
lines. 

Two graphs 

Ti = {C,,Bi,E,,^i), 12 = (£2,^2, ^2, $2) (33) 

are called (topologically) equivalent if there are two maps (pi : Bi —>■ B2 and 02 : 
£1 ^ C2, such that 

^2° 01 = El, (34) 
where o means decomposition of maps, and 



01 : i3i X i3i ^ B2X B2 

01 (^>^) = (35) 

A symmetry of a graph T = (£, B, E, $) is a pair of maps (pi : B ^ B and 02 : £ — >■ 
£, such that 

$ O 02 = 01 O $^ 

Eo0i = E. (36) 

The number of those maps is called the symmetry number of F. 

The set of equivalence classes of connected graphs with E external and L internal 
lines, and with the additional property that every loop has an even number of lines, 
is henceforth denoted by Qe{L), and 

Ge := U ^e{L). (37) 

L>0 

The constraint on the loop is for later convenience, as will become clear below. If 
in addition the graphs are 1-particle irreducible (IPI), i.e. they are still connected 
by removing an arbitrary internal lineQ, we write Qj^^{L) and Q}^^. To any order, 
the latter graphs are considerably less in number. Restricting to even graphs only, 
we define the classes ^^(L), ^^^^ S's'^^^iL) = GW^{L) H Q]^ , correspondingly. 

^We follow standard convention by using the notion IPI (l-"particle" irreducibility). It should 
be kept in mind, however, that the vertex functional r{M,v), eqn. (|^), is not the generating 
functional of the LCE diagrams that are IPI. 
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2.3 Weights and resummations 

The graphical representation of the hnked cluster expansion of susceptibilities is an 
expansion in term of equivalence classes of connected graphs as defined above. For 
every such equivalence class we need exactly one representative. 

Every graph F represents a number, which is called its weight. If F is connected 
and has E external lines, it contributes this weight to the susceptibility xe- For the 
0(N) scalar field models, it is computed along the following lines. Only even graphs 
need to be considered, otherwise the weight vanishes. 

• All the vertices of F are placed at lattice sites. No exclusion principle holds, 
i.e. any number of vertices can be placed at the same site. By assumption, 
the model lives on the lattice A, which satisfies translation invariance. An 
arbitrary selected vertex is located at a fixed lattice site, avoiding a volume 
factor. All the other vertices of the graph are placed arbitrarily. We remark in 
passing that it is only at this point where the geometry of the lattice enters. 
This includes dimension and topology, such as toroidal symmetry for the case 
of a finite temperature lattice. 

• Every internal fine of F is assigned a symmetry label a e {1, . . . , N}. The line 

contributes a factor Vaa{x, y), where x,y ^ A are the lattice sites its endpoints 
are placed at. All the external lines get the same fixed symmetry label attached 
(usually 1 for convenience). 

• A factor for every vertex v & Br- Let 2n be the sum of internal and external 

lines attached to v. Then this factor is given by ^^J^"^|| C2n(Qi, . . . , a2n), where 
ai, . . . , are the symmetry labels of the internal and external lines attached 
to V. 

• The sum is taken over all possible assignments of symmetry labels to lines, 
and over all possible placements of the vertices at lattice sites (except for the 
fixed vertex). 

• Two final factors complete the weight. The first one is given by (S'(r))~^, 
where S{r) is the symmetry number of F. It has the product representation 

S(r)^Sp{r)- n rn(v,w), (38) 

{v,w)eBrxBr 

where the integer number Sp{T) corresponds to the symmetries of F under 
permutations of vertices. The latter product is due to the exchange of multiple 
lines between vertices. The other final factor counts the various enumerations 
of the external lines, 

n.,SrE{v)\- ^^^^ 
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It accounts for the fact that the definition of a graph used here does not give 
identities to the external hues (anonymous factor). 



For the 0(N) models discussed here, summation over all 0(N) symmetry labels 
and the placement of the vertices on the lattice factorizes, i.e. lattice imbedding 
numbers and 0(N) weight factors can be computed independently for every graph. 
Weighted susceptibilities such as the moment fj,2 of ([T^) are described by the same 
diagrams as unweighted ones. The only difference is that they have the appropriate 
weight factors inserted such as g{xi,X2), Xi and X2 being the sites of the external 
vertices. This, of course, must be done before the sum over the lattice imbeddings 
is carried out. 

Placing vertices at lattice sites is subject to strong restrictions by the explicit 
form of v{x,y). Computations have to take this into account explicitly. From now 
on we restrict attention to nearest neighbour interactions only. That is, v{x,y) is of 
the form @). The LCE then is an expansion in powers of 2k. The order of /t a graph 
contributes to is given by its number of lines. Furthermore, the final imbedding of 
the graphs is done onto the hypercubic lattice A. Vertices connected by at least 
one line have to be placed at nearest neighbour lattice sites. A necessary condition 
for a non-vanishing imbedding to exist is that every loop of this graph has an even 
number of line^. It is meaningful to take this into account from the beginning. This 
implies a considerable reduction on the number of diagrams to be considered, as we 
have done with the definition of Qe above. 

Connected graphs and in turn susceptibilities can be explicitly written in terms 
of IPI ones. 

^iPi 

= l-(2«:)(2D)xr' 

+ i2K)i2D^)ixr)' .... 
{I - {2K){2D)xry ' ^ ^ 



(1 - i2K)i2D)xrr ^^^^ 
The IPI susceptibilities admit the series representation 

xT = E«S(A)(2'^)^ (44) 

L>0 

<l(A,a) = E <^Ti\G), (45) 



^We remember that Lq, the possibly finite extent of the lattice in one direction, is always chosen 
to be even. 
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where a^^^{X,G) is the weight of the diagram G that is computed as described 
above, except for the factor {2k,)^ that we have written out exphcitly. Similar 
representations hold with the obvious modifications for e.g. ^l^^. According to the 



identities (^0|f) it is sufficient to restrict attention to IPI graphs only. Actually, this 
simplification implies even more than just a reorganization of diagrams. It turns 
out that critical surfaces are computed from the series expansion of xl^^ with much 
better accuracy than from the series of X2- They are determined by the roots of the 
denominator of (^). 

Further simplification arises by classifying diagrams in terms of even more strin- 
gent properties An essential characterization that we adhere to is 1-vertex 
irreducibility (IVI). A graph F is called IVI if the following condition is satisfied. 
Remove any vertex f of F, together with the external and internal lines attached to 
V. We denote the resulting graph by F^, = B^, Ey, $^). Then, for every w E 
there exists an external vertex u E By and a path from u to w inT^. In other words, 
every connected component of F„ has at least one external line left attached to one 
of its vertices. We write 

Sk{L) = {Te gr\L) I F is IVI } (46) 

and = Sk{L) fl for the set of graphs that are both IPI and IVI. On the 

other hand, renormalized moment diagrams are IPI graphs that have exactly one 
external vertex, 

Qfc(L) = {F G gl^\L) I there is v e Br with Er{v) = k}, (47) 

and Qj^i^L) = Qk{L) fl With both of these notions, susceptibilities are repre- 
sented as 

xT=i: E aT{\,^,G){2n)^, (48) 

where the are computed in the same way as are the a]^^ above with the following 
exception only. The vertex couplings f 2„ are replaced by the renormalized moments 



c _ 

L>0 GeQZiL) 



^2n-^2n=E E ^^(A, G) (2«:)^ (49) 



and the a^^^ are computed as in (^11) before. 

IVI and renormalized moments provide a computational device, which we can 
think of as a kind of product representation of IPI graphs by IVI diagrams and 
renormalized moments. In the final end, after having determined all the weights of 
the graphs, the coefficients are reorganized to give the coefficients a^^i{X). 



3 Algebraic representation of graphs. 
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3.1 Incidence Matrices 



In order to investigate graphs algorithmically and by a computer we need a con- 
venient algebraic representation. Examples are linked list connectivity structures 
or incidence matrices |]TU|, The latter appear to be more useful for numerical 



operations. 

Graphs of the LCE do not have tadpole lines. This allows for the following 
definition. Let F be a graph, Vp the number of its vertices. Then an incidence 
matrix Jr of F is defined as follows. 

• Enumerate the vertices of F, so that 

Br = {v,\i = l,...,Vr}. (50) 
In particular, Vi ^ Vj for i ^ j. 

• For i, j = 1, . . . , Vp, i j, define 

Ir{iJ) = m{vi,Vj), (51) 
where m{vi,Vj) is the number of lines between Vi and vj. 

• The diagonal elements are given by the number of external lines attached to 
the vertices, 

Ir{i,i) = E{v,), 2 = l,...,Fr. (52) 

The incidence matrix is a symmetric Vp x Vp matrix, Jp G M(Vr x Vp). Most 
algorithms need to deal with e.g. the right upper triangle only. 

The construction above does not provide a unique definition of an incidence 
matrix of a graph. Any other enumeration of its vertices correspond to a non-trivial 
permutation vr E Uy^ and leads to the incidence matrix 

J?(z,j) = /r(7r(^),7r(j)), (53) 

i.e. Jr is unique only modulo simultaneous permutations of its rows and columns. 
Those permutations that leave Jr unchanged constitute a finite subgroup of Ilvr. 
The number of elements of this group is equal to S'p(F). Two graphs Fi and F2 are 
topologically equivalent if and only if they have the same number of vertices V, say, 
and their respective incidence matrices (with respect to an arbitrary initial labelling 
of the vertices) satisfy 

Ir. = It,, (54) 

for some tt G fly. 

During the process of graph generation discussed in the next section, many equiv- 
alent graphs will be constructed, but for every equivalence class we have to keep only 
one graph. To decide whether two graph are equivalent according to (^^ will need 
factorials of permutations to be carried out. This is too expensive for high order 
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graphs. What is needed is a unique representation of graphs by incidence matrices. 
Furthermore, the computation of those incidence matrices should be as efficient as 
possible. Definitely huge factorials have to be avoided. 

The representation of graphs by incidence matrices can be made unique if a total 
order relation on M{V x V) is chosen for every V, and in turn the maximum over 
all simultaneous permutations of rows and columns is taken. Then, for an arbitrary 
initial labelling, we get a canonical representation by 

/r := max Jp, (55) 

and the maximum is with respect to the total order relation. A convenient choice 
is as follows. Let A, B e M{V x V). Then 

A>B& there is i,j e with > and 

A{k, I) = B{k, I) for a\\k = i,l < j and for all k < i. ^ ^ 

Still there are huge factorials involved to construct the canonical form of incidence 

matrices. The breaktie solution here is to introduce an order relations for the vertices 

of a graph first. Topological properties of the graph are taken into accont here, local 

and less local ones, as we will discuss below. Ordering cannot be complete, at 

least those vertices need to stay unordered relative to each other whose exchange 

correspond to a symmetry of the graph. Once the order relation has been chosen, 

a canonical incidence matrix can be defined by subsequently allowing in (|55|) only 

those vertices to be exchanged that stay unordered relative to each other. 

For any graph F let O be an order relation on its set of vertices Br- For any pair 

o o o o o 

of vertices v.w E Br, v w, either v>w, w>v, or t>~ti; if neither v>w nor w>v. The 

o o o o 

enumeration of the vertices then is any one such that Vi>~V2>~V2,>~ ■ ■ ■ >~vy^, where 
o o o 

>- means > or ~, leading to some Jr. The canonical incidence matrix of F then is 
defined by 

I? := max /f , (57) 

with Ily^ the subset of IIvp whose elements correspond to the exchange of those 
vertices only that stay unordered relative to each other with respect to O. Further- 
more, the number of permutations tt with = /f is equal to the (partial) symmetry 
factor Sp(V). 

This definition of an incidence matrix satisfies all the above criteria. The partic- 
ular form depends on the order relation O. The expense of ordering vertices is with 
some power of the number of vertices, much cheaper than factorials. 

The heart of the problem is to find an appropriate order relation O on the vertices 
of a graph. The number of remaining permutations of (|57|) should be as small as 
possible. Typically it is still some factorial, but smaller than the original one. To 
some order of the LCE we succeed with a chosen vertex ordering, i.e all the graphs 
of a particular correlation function can be generated with acceptable costs. With 
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increasing number of lines, however, the problem of large factorials reappears, and 
one is forced to extend the order relation. The cardinality of Uy is bounded from 
below by the partial symmetry number Sp{r). There is much space left towards 
this ideal case. 

For later convenience we introduce the following definition. Let Oi, i = 1, . . . ,n, 

be order relations on the set of vertices Br of a graph F. Then we define an order 

relation O =< Oi, . . . , 0„ > on i3r as follows. For v,w E Br, v is called larger than 

o 

w with respect to O, v>w, if there is z G {1, . . . , n} such that 



v>w (58) 



and for all j < i 



(59) 



If for alH G {1, . . . , n} f ~w, v and w stay unordered with respect to f We 
notice that the order relation O depends on the order of the sequence Oi, . . . , 0„. 

In connection with the LCE, the first order relation on the vertices is due to 
Liischer and Weisz 0. It takes into account all local properties of the graph in 
an optimal way. That is, the order number assigned to a vertex v is based on the 
connections to its neighbours N'{v). Let Vi{v) denote the number of vertices of a 
graph F that have precisely i lines in common with the vertex v. Then 

=< Oa,0/3 >, (60) 

where for all f , G i3r, 

v>w ^ E{v) > E{w), (61) 

Of3 

v>w <^=^ there is i > 1 with z/i(f) > Uii^w) and 

i^j{w) for all j > i. (62) 



By means of this order relation the graphs of the two- and four-point correlations 
can be generated up to and including the 14th order graphs. This order appears as 
some natural limit, due to the kind of slowing down mentioned above. As we have 
discussed in Sect. 1, in order to obtain most accurate quantitative information on 
critical behaviour from the LCE series, we want to compute them up to the 18th 
order. In particular we expect to be forced to do so on finite temperature lattices. 

The problem that arises becomes clear already by considering the 18th order 
graph F of Fig.|I|. The original 17! number of permutations to be carried out in 
order to construct the incidence matrix according to ( p5D is reduced after Olw to 
15!, which still is a large number. This is because the graph has a large number of 
2-vertices. On the other hand, there are obvious subsets of 2-vertices whose permu- 
tations cannot be part of a symmetry. Those vertices should be further distinguished 
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Figure 1: Motivation for extended vertex ordering. The graph T shown has many 
2- vertices. They stay unordered relative to each other by the order relation O^w 



by the order relation on Br- For instance, as an immediate criterion we can use the 
distance of such a vertex along a chain of 2-vertices towards the next non-2-vertex. 
Before we discuss the extension of the order relation Olw, we notice that for 

every pair of vertices v,w & Br with v G B^lntj w ^ B^lnt^ we have either v >' w 

or w > V. That is, Olw separates internal 2-vertices from those that are not. 
Extension of Olw can thus be done for internal 2-vertices and their complement 
separately. 

Henceforth we need the following definition. Let Vi, V2 C Br be complementary 
subsets, Vi n V2 = 0, Vi U V2 = Br, and let Oi and O2 be order relations defined on 
Vi and V2, respectively. Furthermore, let Oq be an order relation on Br such that 

Co Oo 

for all e Vi, ^2 G V2 either vi>V2 or V2>vi. Then we define an order relation 
O =< C»o, O2] > as follows. For all v,w e Br, 



o 

v>w 



^0 • r M -11 

v>w , it f G w G 1/2 or vice versa, 

v>w 



, ifv,weVi, (63) 
v>w , if v,w E V2, 



00 o 
where 0[ =< Oq, Oi >, =< ^o? O2 >■ If neither v>w nor w>v, we write 

This can happen only if v, w both belong to Vi or both to V2. 

3.2 Extended vertex ordering for non-ring graphs. 

In the following we restrict attention to connected graphs only. We have to distin- 
guish between ring graphs and their complement, non-ring graphs. By definition, a 
connected graph F is a ring graph if F has at least three vertices, and Br 7^ Br int- 
That is, F has at least one vertex that is not an internal 2-vertex. We discuss ex- 
tended vertex ordering of non-ring graphs first, the modifications for ring graphs are 
described afterwards. Let us suppose in this section that F is not a ring graph. 



We define the order relation O used with (57) to define the canonical incidence 
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matrix by 

O =<OLW,[Oni.Oint]>, (64) 

where Oint and Oni are convenient order relations on the set of internal 2-vertices 
and its complement, respectively. As we shall see they can be chosen in such a way 
that the vertex ordering is considerably enhanced. We are now going to define Oni 
and Oint- 

A path p = {wi,W2, . . . , Wn) is called a 2-chain of wi if W2, . . . , G i3|?]„t are 

mutually distinct. If the last vertex is not an internal 2-vertex, Wn ^ ^rint^ P is 
called maximal 2-chain of Wi. In this case we assign a tripel of two integers and a 
vertex to p by 



M = n-1, (65) 

n-l 

Kp) = '^ni{wi,Wi+i), and (66) 
1=1 

e(p) = Wn- (67) 

For any v ^ B^^ there are precisely two distinct maximal 2-chains. Let us denote 
them by ci{v) and C2{v). They intersect in v and possibly the last vertex only. For 
later use we set 

/i(t>) = max /i(Q). (68) 

1=1,2 
(2) 

Let Onon be an order relation on Br \ Bf^jnt- induces an order on the maximal 
2-chains {ci(f ), C2(f )} as follows. For i,j e {1,2}, i ^ j, 

fx{ci) < fx{cj) , or 

Ci{v) > Cj{v) ^ Kci) = Kcj) and A(q) > X{cj) , or (gg) 

yu(cj) = /i(cj) and A(cj) = X{cj) and e(cj) >"e(cj). 

If unique, we choose " short and long chains" by 

Cs{v) = max Cj(f), (70) 

i=l,2 

ci{v) = mill Ci{v). (71) 

Roughly speaking, Cs{v) gives a stronger binding of v to e(cs) than Ci(i;) to e{ci) 
(shorter chain, more lines). If neither Ci{v) > C2{v) nor C2(f) > Ci{v), we choose 
any one of the maximal 2-chains as Cs{v) and the other one as ci{v). This freedom 
is allowed because all that will be used is the order relation (^^. We remark in 
passing that e(ci) and e(c2) can be identical. 

Finally, Onon induces through the definition (|69D an order relation Oint on B^l^^ 
by Oint =< Oii,Oi2 >, with 

v>w Cs{v)>Cs{w), (72) 

t;>w; ci{v) > ci{w), (73) 
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for all internal 2-vertices v,w. This is defined once Onan is fixed. One possibility 
is to define Onon as equal to Olw- But it appears that only part of the rela- 
tions lead to considerable enhancement of ordering, other ones in very exceptional 
cases only. Furthermore, in addition to supplementing the ordering of internal 2- 
vertices, we also extend the order relations on the other vertices. Possible criteria 
are properties of attached 2-chains. Let us call such a relation Oni- Then we let 
Onon =< Oa, 0^,0s, Oni >, whcrc Oa is defined by (l6l|) , and 

v>w N{v)>N{w), (74) 

v>w <^==^ max k > max fc, (75) 

for all v,w e Br\ B^jnt- 

AU that is left to be done is to define Oni- In the first place Oni supplements Olw 
for vertices that are not internal 2-vertices. We take into account properties such 
as the lengths of attached 2-chains, i.e. their number of vertices, and the number 
of lines thereof. In addition to these direct enhancements we iterate the ordering. 
That means, having fixed order relations in a particular way, we can impose even 
more ordering by taking into account order criteria of neighbours. To this end, we 
define the following numbers. For every v E BrXB^ 



inti 



d,{v) = E{w), (76) 

w£j\f{v) 

d,iv) = E Niw), (77) 

T,xnt 

d4iv) = /^(^)' (79) 

T,%nt 

d^{v) = max fi{w). (80) 

T,%nt 

Then the final definition of Oni is as follows. 

Oni ^ Oni^lj - - - j Oni,5 >, (81) 

where for all /c = 1, . . . , 5 and for all v,w E Br\ B^lnt, 

V >''w dk{v) > dkiw). (82) 

After all these technical details let us now give some examples. To compute the 
canonical incidence matrix of a graph F under consideration, the vertices are 
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ordered first according to the order relation (0). In turn if is determined by (|57D. 
For the graph of Fig. H above there are (3!)^ permutations left, orders of magnitude 
less than 15! or even 17!. Two further examples are shown in Fig. 0. Those diagrams 
belong to the intermediate graph classes out of which the "physical" graphs such as 
iSfc and Qk are constructed, as will be discussed in the next section. 

Figure 2: Two examples for the efficiency of extended and iterated preordering. 
The bottom lines below the graphs show the number of permutations left after 
no preordering at all, preordering according to Liischer and Weisz, and extended 
iterated preordering. 




16! 4!10! ^ 2^ 14! ^4- 10! ^2^ 



3.3 Ring graphs. 

A special treatment is necessary for those graphs that have internal 2-vertices only. 
Let r be such a ring graph, so that Br = B^lnt- The distinction between internal 
2-vertices and those that are not for non-ring graphs is replaced now by identifying 
so-called symmetric 2-vertices. Such a vertex is defined by the property that with 
both of its neighbours it has the same number of lines in common. More precisely, 
a vertex f G -Br is called a symmetric 2- vertex if there is an integer G N such that 

uj,{v) = 2. (83) 

The number k is denoted by c{v) (There is at most one such k for every vertex. For 
/ 7^ k, viiv) = 0). The set of symmetric 2-vertices is denoted by Br^symm- A ring 
graph r is called exact if every vertex f G -Br is a symmetric 2-vertex. 

The definitions and constructions of 2-chains, order relations on them and sub- 
sequently on internal 2-vertices as discussed in the last section for non-ring graphs 
can now be translated to ring graphs. The essential step is to replace the notion 
of internal 2-vertices by symmetric 2-vertices, and similarly for their complements. 
Some of the notions used to define the order relations above have to be adjusted 
correspondingly, other ones become dummy. The details are very similar as in the 
last section. For completeness they are given in the Appendix for non-exact ring 
graphs. 
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We are left to define the canonical incidence matrix for exact ring graphs. This 
is done explicitly. Let F be an exact ring graph with L lines and V vertices. Then 
y > 3, and v = L/V a. positive integer. The canonical incidence matrix Jr is 
defined as the symmetric V x V matrix with the only non-vanishing entries of its 
right upper triangle being given as follows. 

7r(fc,fc + 2) = u , hr k = l,...,V -2, (84) 
7r(l,2) =7r(\/-l,y) = u. (85) 

In particular, the diagonal entries are zero. The symmetry number is given by 
^(r) = 2V-{iy\)^. 



4 Generation of graphs and weight computation. 

In this section we present and discuss efficient algorithms to generate all the equiv- 
alence classes of graphs we are looking for to compute correlation functions and 
susceptibilities. The central idea is to generate them from base classes of diagrams 
that are simpler in structure and less in number. These classes should be universal 
in the sense that the various correlations are generated from the same base classes. 
The latter have to be provided just once. The actual order by order construction of 
graphs on the number of lines is essentially done in the simplest class only. 



4.1 Inheritance tree of graph classes 

In Sect.0 we have outlined the way connected correlations and susceptibilities can 
be built up from IPI ones. The latter allow for some kind of product representation 
by IVI graphs Sk and renormalized moments Qk, (^^-^). Only those graphs have 



to be constructed explicitly. The most apparent difference between graphs of Sk 
and Qk is the way they have external lines attached. Every F G iS^ has at least two 
external vertices, whereas all external lines of any F G enter the same vertex. 

We generate those diagrams from general IPI ones with internal lines only by 
attaching external lines in the various different ways. This is possible because the 
IPI property is independent of external lines. Even further simplification arise if 
we omit multiple lines between vertices. Any IPI graph with multiple lines can be 
constructed from a connected graph with single lines only, by multiplying internal 
lines appropriately. Actually, care is needed here because we have to relax the IPI 
property for those simplest graphs. 

Fortunately, we can work around this problem. We allow for double lines under 
special circumstances, still having an acceptibly small base class. This design frees 
us from having to check IPI properties for more advanced classes during the con- 
struction process. Also, graphs do not need to be glued together. We are now going 
to state our design of graph classes and the sequence of construction mechanism. 
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Although the description here is in terms of graphs themselves, the implemen- 
tation on the computer is done in terms of incidence matrices and is more or less 
straight forward. Let 

r = (£r,i3r,^r,$r), 
be a diagram and C £p a subset of the lines of F. Then we set 

T\M = (Cr,l3r,Er,^r), (86) 
where Cr = '^rX-M. and $r is the incidence relation $r restricted to ^r, 

$r = $r|zr • 

r\A4 is the graph F with the lines of Ai omitted. With this notion, the base classes 
Vi (L) for the various numbers of lines L are defined as follows. 

Vi (L) is the set of graphs F that have L internal lines, no external lines, and are 
IPI, i.e. 

F e C'(L). (87) 

Furthermore, for all pairs of vertices v,w ^ Br, either there is at most one line 
between v and w, m{v,w) < 1, or m{v,w) = 2 and F \ J^{v,w) is not connected. 
Here, J^{v,w) is the set of lines between v and w, as defined by (^91) . 

Figure 3: The structure of the graph class Vi. The hatched bubbles represent the 
reduced components, which are connected by double lines. 



Thus, removing any one of a double line between two vertices destroys the IPI 
property, removing both disconnects the graph. Fig. ^ shows the typical structure of 
a graph F G Vi{L). Removing all double lines, F decomposes into various connected 
components. The latter are IPI, and any two vertices have at most one line in 
common. Those parts are called reduced components. Replacing each reduced 
component by a single vertex, and every double line of F by a single line, the 
resulting graph is a tree, i.e. it does not have any loop. 

From the Vi we derive the next more involved and larger classes V2{L) which still 
lacks external lines, but have an arbitrary number of lines between their vertices, 

V,{L)=gr. (88) 
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All the graphs of V2{L) are obtained by taking every graph of Vi{L') and multiplying 
its lines in all possible ways, so that in total L — L' lines are added. This has to 
be done for every 2 < L' < L whenever L > 0. In this way, for every equivalence 
class of 7^2 (-^) there will be at least one graph constructed, i.e. 7^2 (-^) is covered. 
However, many graphs generated in this way are equivalent, and we have to keep 
track that for every equivalence class only one representant is kept. Fortunately, we 
don't have to compare a newly generated graph to all prior constructed diagrams. 
This is due to the fact that two graphs of V2{L) can be equivalent only if they are 
both generated from the same graph of some Vi{L') according to the above rules. 

The last step of the graph generation sequence is to generate the IVI graphs 
iSfc(L) and the renormalized moments Qk{L) out of the graphs of V2{L). This is 
done by fetching all graphs of V2{L), and attaching k external lines in different 
ways. Attaching all of them to just one vertex generates Qk{L). 

Figure 4: Inheritance tree of graph classes Vi V2 ^ S and Q. The right hand 
side of the figure shows an example construction. Actually, the graphs of Q^^^ are 
no more constructed explicitly. 




The graphs of Sk{L) have at least two external vertices once L > 0. Furthermore, 
they have to be IVI, which restricts the number of possibilities to attach the external 



lines. To check this property needs one recursive depths first search ||Tl| within the 
augmented graph. It is obtained by adding one vertex and attaching all external 
lines to it f^, |12| . This is a very cheap algorithm and can be done before the incidence 
matrix of the graph is brought to its canonical form. A further simplification arises 
for the 0(N) scalar models due to the fact that we can restrict Sk and Qk to even 
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graphs only. 

Again, only graphs that are mutually non-equivalent have to be kept. Similar as 
before, two graphs can be equivalent only if they originate from the same graph of 

The generation sequence of the classes is scetched in Fig. |[ 

It remains to construct the base classes Vi{L). This is done inductively on L. 

We closely follow the proposal of |Q, adjusted to the present case. The lowest order 

graphs are given explicitly, Vi{L) for L < 4, as in Fig. ||. 

Figure 5: The graphs of Pi{L) for all L < 4. 
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Now suppose that L > 5. Then Vi{L) is constructed from Vi{L') with L' = 
L — 1, L — 2 by three rules, 

Vi{L-l) ^ Vi{L) , Pi(L-2) "^^^ Vi{L). (89) 

They apply to every graph of the respective class and operate separately on every 
pair of vertices that satisfy appropriate conditions. To state these rules we have to 
provide some definitions first. 

Let r G Vi{L) for some L. Then to every vertex v E B-p we assign 

• a parity number 77(f) G { — 1, +1} in such a way that neighboured vertices have 
opposite parity. This is unique up to a global sign for F, due to the fact that 
every loop of F has an even number of lines. 

• a reduced component index r{v) G {1,2, ... ,n}, where n denotes the number 
of reduced components of F. This is done in such a way that vertices that 
belong to the same reduced component get the same number r attached. Those 
numbers are unique up to permutations of {1, . . . , n}. 

With these notions, the rules operate as follows. 

Ri on F G Vi{L — 1) and for Vi,V2 G Br with 
Vivi) 7^ r]{v2), m{vi,V2) = 0, r(fi) = r{v2), 
insert a line between vi and V2, 



^ ' . (90) 

Vi V2 Vi V2 



In particular, Vi 7^ ^2- 
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on r G — 2) and for Vi,V2 G By with 

r]{vi) = r]{v2) (so that m{vi,V2) = 0), r{vi) = r{v2), 
insert a vertex w and 2 single hnes connecting w to Vi and to f2, 



• ^ • • • • (91) 

Vl V2 Vl W V2 

There is no hne and there will be no line between vi and V2- 

i?3 on r G Vi{L — 2) and for ^1,^2 G Br with 

■M{vi, V2) 7^ 0, and at least one vertex has more than 2 lines attached. 

Let I G M.{vi,V2)- This line I is replaced by 2 vertices 1^1,^2 and 3 new lines 

forming a 2-chain (f 1, ifi, W2, ^2) of single lines, 



■ (92) 

Vl V2 Vl Wi W2 V2 

Some remarks are in order. For the first two rules, both vertices must belong to 
the same reduced component, otherwise the class Vi would be left. In the 2. case, 
the two vertices Vi and f 2 can be equal. In this case a double line is generated, and 
a new reduced component appears, consisting of a single vertex only. 

Figure 6: To a long 2-chain, rule R3 need to be applied only once. 



The 3. rule applies also if the line / belongs to a double hne. Obviously, we need 
to operate only on one of those double lines. Two reduced components melt to a 
single one in this case. Furthermore, for a 2-chain as in Fig. it is sufficient to 
operate on the left or right end. 

Again we have to check that every newly generated graph provides a new equiv- 
alence class. In contrast to the construction of V2, S and Q discussed above, we 
have to compare to all the graphs generated yet (with same number of lines and 
vertices, of course). But this number is considerably smaller than that of the more 
involved classes. 

Definition and construction mechanism of the intermediate classes is not unique, 
of course. Our choice and design of the classes differ from those of Liischer and 
Weisz [|]. We allow graphs of the base class Vi not to be twice connected. This 
increases its number of graphs, but it is still acceptibly small. The major advantage 
of our approach is that we do not have to glue graphs in any way. The costs increase 
with minor power with respect to the number of graphs. 
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4.2 Weight computation. 

Every graph represents a real number it contributes to a susceptibility, its weight. 
For the 0(N) scalar models the most laborious parts of it are the computation of 
both the lattice imbedding number and the 0(N) factors. The other contributing 
factors are more easily obtained. In particular, the topological symmetry number is 
computed in the course of the graph generation. 

The lattice imbedding number of a connected graph F is given by the number of 
possibilities all its vertices v & Br can be placed at lattice sites x{v) consistent with 
the topology of F. That is, vertices connected by a line have to be placed at nearest 
neighbour lattice sites. One vertex can be fixed, saving the extensive volume factor. 
It is the nature of the LCE that no exclusion principle applies, i.e. different vertices 
can be placed onto the same lattice site. 

For the explicit computation of the lattice imbedding number we use the tech- 
niques proposed in [|| appropriately generalized. The central idea is to place only 
those vertices v explicitly that are not internal 2-vertices, v ^ B-^int- They are 
connected by maximal 2-chains. The latter are considered as representing random 
walks of fixed length. Once all v ^ B^lnt ^-^^ P^t onto the lattice, every 2-chain 
contributes a factor given by the number of random walks of its length from one 
end point to the other. Closed analytic expressions can be derived for them. In 

('2) 

turn, summation is done over all placements of the v ^ B-p i^^^. The shortest 2-chain 
between two vertices determines the maximal distance of the lattice sites they can 
be placed at. On the Lq x oo^^^ lattice Ai, the distance is defined by 

D-l 

{x,y) = \\xQ-yo\\Lo+^\xi-yi\ iorx,yeAi, (93) 

i=l 

where ^ 

IIxoIIlo = inf |xo + nLo| < (94) 

raGZ Z 

(2) 

To avoid unefficient trial and error, placing the vertices f ^ i3p on the lattice 
is done in fixed order determined by the so-called imbedding preordering. This 
enumeration is defined as follows. The only vertex chosen to be placed at some fixed 
lattice site is attached the smallest number 1. All the vertices v ^ -Sf^Lt attached to 
it by a 2-chain, say m in number, get the numbers 2, . . . , m + 1 attached, in arbitrary 
order. In this way we proceed iteratively. Having enumerated a subset of vertices, 
enumeration is continued with those vertices v ^ S^lnt that as yet have not been 
attached an imbedding preorder number, and that are connected to the subset by a 
2-chain. This is done until all v ^ have got a number. One breadth first search 
is necessary to obtain this enumeration. The actual placement onto the lattice is 
done in the order of increasing imbedding preorder number. 

Connectivities have to be respected. This is best done by the imbedding pre- 
ordered incidence matrix Jp . Its {i,j)th entry is given by the length of the 

(2) 

shortest 2-chain between the vertices Vi and Vj, both Vi,Vj ^ -Spj^^, with respect to 
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the enumeration just defined. In contrast to the graph construction, we do not have 
to worry about uniqueness. 

Imbedding the graph onto the lattice now proceeds by sweeping over the right 
upper triangle of /p from left to right, top to bottom. Entering a non-zero entry 
Jp (/e, /) implies either that vi has to be placed onto the lattice, or if placed already, 
verify that 

ixivk),xivi))<l^ikj). (95) 
If satisfied, we continue. Otherwise, we re-place vi within the appropriate cube 

ixivi),xM)<I^ilJ'), (96) 

where I' is the smalUest imbedding preorder number such that vi and vi' are con- 
nected by a 2-chain. At the same time, all Vm with m > I are put off the lattice. 
A succesful sweep implies that a particular imbedding is found. All the random 
walk numbers of the 2-chains are multiplied together. Imbedding search continues 
by moving the latest vertex that can be placed to a position not yet taken within 
its cube, and proceed as before. If it is the first (fixed) vertex, we are done. 

The discussion above concerns the computation of the lattice imbedding number 
of a single connected graph. Actually, for classes of graphs we don't need to compute 
all those numbers separately. This is because imbedding numbers do not depend 
on multiple lines between vertices. Only the nearest neighbour constraint matters. 
Replacing multiple lines by single lines, many graphs become topologically equivalent 
and thus have identical imbedding number. To make this working needs canonical 
ordering and buffering. Further improvement is obtained if no explicit weight factors 
such as g{x, y), (0f), are involved. In this case we do not need to distinguish between 
internal and external vertices. 

The way we have done the computation of the 0(N) factors of a graph F, i.e. the 
summation over all 0(N) symmetry labels attached to the internal lines of F, pro- 
ceeds as follows. All the external lines are connected in pairs by (possibly trivial) 
paths of lines through the graph. Once all the external lines are paired in a particu- 
lar way, all the remaining internal lines have to be combined to loops in all possible 
ways. Every succesfull covering of the lines yield a N'^^ , where M is the number of 
loops. The sum over all those configurations yields the 0(N) factor of F. 

Let us make the notions more precise. First of all, to avoid overcounting, we 
enumerate the external lines. Let F = (£p, i3r, -Er, ^r) be a graph, L and E = 
J2veBr its number of internal and external lines, respectively. An enumeration 

of the external lines of F is a map 

r] :{!,..., E}^Br,e.t (97) 
such that for every v G Br,ext, the set 

{te{l,...,E}\r]{i)=v} (98) 
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has precisely Er{v) elements. In practice, enumeration is done in conformity with 
canonical ordering. That is, for i < j, T]{i) < ri{j) or rj{i) ~ rjij)- A pairing of the 
E external lines is a set oi E/2 pairs 

{(ei,e|+J,...,(e|,eE)}, (99) 

with ei, . . . , e^; G {1, . . . , E} mutually distinct. The set of all pairings is denoted by 
Vv{E). 

A path of lines of F is a sequence 

P= {Wi, h,W2, h,---, Wn-l, ln-l,Wn) (100) 

with 

wi,...,Wn e Br , 

h, . . . , /„_i G Cr , mutually distinct, (101) 
Wi, Wi+i e $r(/i), , Wi 7^ Wi+i for all i=l,. . . ,n-l. 

If Wi = Wn, p is called a loop of lines. The set of loops of T is denoted by Or- Two 
path of lines are called disjoint if they don't have any line in common. is the set 
whose elements are the sets of mutually disjoint paths. Furthermore, the ordered 
tripel (ei,p, 62) with 61,62 G {1,...,E}, 61 7^ 62 and r]{ei) = Wi, 77(62) = w„, is 
called an external path of lines. This does not exclude p itself to be a loop. For, 
77(61) = 77(62) is not excluded. £r is the set of all external paths of F. Finally, a set 
of paths {pi, . . . ,Pm}, m < L, is said to cover F if every internal line of F belongs 
to exactly one path. 

With these notions, the 0(N) symmetry number of F is explicitly given and 
determined by 

c{r)= E E (102) 

Ps-Pp{£;) , }e-Dp 0<M<L gi ,...,<jM60r 

T {Pl,...,p^ ,(Jl,...,l3J^J-} covers r 

(f^i.Pi.e E , .)S£r ^ill (si'E E I .)eP,i=l,....S- T 
-7 + ' T+' 

The actual recursive implementation of this formula from left to right is straight 
forward^. Paths and loops are constructed by moving systematically through the 
graph. Many simplifications can be taken into account for multiple lines between 
vertices or multiple external lines at the same vertex. 

As is the case for the computation of the lattice imbedding numbers, here also 
we don't have to compute the 0(N) factor for every single graph of a class. The 
first obvious observation to be made is that the 0(N) number does not depend on 
the length of 2-chains. All the graphs that become topologically equivalent after 
shrinking all 2-chains between non 2- vertices to a single line have identical 0(N) 

■^This is one of the cases where the description of an implementation is more comphcated than 
the implementation itself. 
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Figure 7: Examples of graph reduction in the course of the computation of the 
0(N) symmetry numbers. In the first case, both /i and v are odd numbers, and 
both are even in the second example. 




V 



g{N,i2,u) K 



number. Other reductions are applied by explicitly simplifying parts of the diagram 
that allow for closed computation. Examples are given by Fig. |^. Once a graph F 
has been reduced to a graph F, its 0(N) factor C(F) is the product of 

• the product of all / and g obtained during reduction, denoted by fr{N), and 

• the 0(N) factor C(T) of the reduced graph F. 

Because many graph by reduction result into the same reduced graph F, it is 
convenient to buffer F and C(F). When a later graph F' of the class satisfies F' = F, 
then 

CiT') = friN)-CiT), (103) 
i.e. C(F) does not need to be calculated a new. 

5 Results. 

We summarize the methods and the results obtained on susceptibility series for the 
0(N) scalar field models on the lattice, as described by the partition function (|l|)- 
(^. The lattice is assumed to be of hypercubic geometry, with periodic boundary 
conditions imposed for those directions in which the extent is finite. Fields at 
different lattice sites are coupled by nearest neighbour pair interactions of strength 
2k. 

Susceptibilities are given as power series with respect to the coupling strength 
2k, with coefficients depending on the bare coupling constants A. These series have 
a non- vanishing radius of convergence. The coefficients are represented as sums over 
equivalence classes of connected graphs with appropriate weights. The number of 
internal lines of a diagram is equal to the order of 2k it contributes to. 
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Table 1: The number of ineqiiivalent graphs as yet generated of the various graph 
classes defined in the last sections. L denotes the number of lines. The renormalized 
moments have the same number of graphs independent of the integer n. All 
the graphs are IPI. Both the renormalized moments and the IVI graphs »S are even, 
i.e. every vertex has an even number of lines attached. 



L 


Pi(L) 


V2{L) 


Q7{L) 











1 


1 


1 


1 


1 


1 


1 




















2 


1 


1 


1 





1 


1 


3 





1 





1 


1 


2 


4 


2 


3 


4 





4 


6 


5 





3 





2 


4 


11 


6 


5 


11 


15 


3 


20 


46 


7 


1 


16 





8 


27 


91 


8 


15 


53 


79 


9 


117 


349 


9 


7 


112 





40 


214 


837 


10 


57 


354 


439 


68 


815 


3140 


11 


48 


953 





247 


1830 


8401 


12 


278 


3160 


2877 


470 


6721 


31187 


13 


379 


9909 





1779 


17028 


90599 


14 


1647 


34457 


20507 


3937 


61653 


336582 


15 


3328 


119921 





14801 


170923 


1042392 


16 


12321 


439552 


161459 


35509 


621191 


3895341 


17 


31869 


1638878 





135988 


1834324 




18 


111493 


6312209 


1376794 


350614 


5548427 




19 


337817 


24810845 





1361878 


20967387 





With increasing order, the construction and computation of all the connected 
graphs of a correlation function or susceptibility becomes managable only if strict 
bookkepping is done. Both the number of graphs and their complexity are consid- 
erably reduced by applying partial resummations. The first essential step in this 
direction is the representation of connected correlations in terms of one-particle ir- 
reducible (IPl) ones. The latter are sums of IPI graphs only and have similar series 
representations, 

xT = E <IW (2«^)^ (104) 

L>0 
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and similar for weighted correlations such as 

/^^ = E/^sw(2'^)^ (105) 

L>0 

cf. (p!2D , (^TD ■ For instance, 2-point susceptibilities are related by 

cf. (^0|f). A further useful reorganization is some kind of product representation of 
weights in terms of graphs that are in addition one-vertex irreducible, and so-called 
renormalized moments, the classes S and Q. It is sufficient to construct the latter 
ones explicitly. The expensive computation of the weights of graphs must be done 
for those classes only. 

In order to manage diagrams numerically, every equivalence class of graphs is 
represented by an incidence matrix. This representation is made unique by defining a 
total order relation on the set of symmetric matrices, and than invoke an appropriate 
maximization procedure over all simultaneous permutations of rows and columns. 
Those permutations correspond to different enumeration of vertices. To avoid huge 
factorials of permutations to be carried out, first of all order relations on vertices 
are introduced, leading to a preordered incidence matrix. The canonical incidence 
matrix then is defined by maximizing only over those permutations which exchange 
vertices that are left unordered relative to each other. 

It is at the heart of the construction mechanism to make this representation as 
efficient as possible. We have to introduce an extraordinary far-going ordering on 
vertices to obtain all the graphs contributing to the susceptibility series up to 18th 
order and beyond. The numbers of the graphs as yet generated are listed in Tab. |1|. 
For the 0(N) scalar models, only those graphs have non- vanishing weights that have 
even number of lines at every vertex. In turn the number of external lines is even 
as well. 

Storage of that many graphs or incidence matrices is as binary words, using a 
kind of HufFmann tree coding mechanism. It is done in such a way that all graphs 
with same number of internal, external lines and vertices have identical bit lengths. 
This is most convenient for the many cases graphs have to be compared. What is 
compared actually are their binary words. 

The computation of the weights of the graphs is done as outlined in the last 
section. The laborious part here consists of the determination of lattice imbedding 
numbers and 0(N) symmetry factors. In total, two orders increase the number of 
graphs of the classes considered by a factor of about 10. Increasing complexity 
implies a total cost factor of about 30-40 to compute two orders more. 

Two final problems have to be managed. Vertices contribute to the weights in 
dependence on the bare coupling constants. In principle, summation of millions 
of graphs must be done each time a coupling constant changes. Even worse, the 
weights typically are of comparable size, so we have to be aware of severe roundoff 
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Table 2: Examples for the number of vertex structures of various graph classes. 



L 


QT{L) 


QTiL) 




sr{L) 







1 


1 


1 


1 


1 


1 

















2 


1 


1 





1 


1 


3 








1 


1 


2 


4 


4 


4 





3 


4 


5 








2 


4 


8 


6 


10 


11 


3 


9 


14 


7 








6 


10 


20 


8 


21 


27 


7 


19 


29 


9 








11 


24 


42 


10 


42 


50 


19 


39 


60 


11 








28 


48 


81 


12 


78 


102 


39 


74 


110 


13 








49 


92 


147 


14 


139 


173 


78 


135 


198 


15 








102 


169 


259 


16 


242 


312 


146 


238 


338 


17 








178 


297 




18 


408 


505 


264 


404 





errors. The total weight of a graph consists of two factors. The first is a non- 
negative rational number, itself a product of lattice imbedding number and 0(N) 
and topological symmetry numbers. The second one comes from the product over 
the vertices. Every vertex contributes a finite dimensional integral that in most 
cases can be done numerically only. These integrals alternate in sign with respect 
to half the number of lines entering the vertex. 

We solve the problem of both having to do a huge sum each time the couplings 
change and the accumulation of roundoff errors by introducing so-called vertex struc- 
tures. A vertex structure maps a graph onto a multitupel of non-negative integers, 

V : ({0,1,2,...})", (107) 

r (//i,//2, ...). 

Here, n is an upper bound on the number of the sum of internal and external lines 
of a graph under consideration, and /li denotes the number of vertices of the graph 
r that have precisely i lines entering, i.e. of 

{v eBr\ liv) = i}. (108) 
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The multitupel (/ii,/i2, . . .) is called the vertex structure of T. For the 0(N) case, 
only /ij with even i can be different from zero. 

Table 3: Coefficients of susceptibility series for the 0(4) model at Ai = cxd and 
A2 = 0, on the 4 x 00^ lattice. 



L 


1 -1 TST 1 


1 i TST 1 


1 -1 TST 1 


1 -1 TST 1 

IPl 





0.25 





-0.0625 


0.078125 


1 














2 


-0.125 





0.1875 


-0.6640625000 


3 


0.01041666667 


0.0078125 


-0.04166666667 


0.2669270833 


4 


0.01627604167 





-0.1513671875 


1.236368815 


5 


0.001139322917 


0.0004882812500 


0.01839192708 


-0.3427937826 


6 


-0.03774685330 


0.002115885417 


0.1999131944 


-2.336894565 


7 


0.008895534939 


0.006170654297 


-0.08401557075 


1.321710798 


8 


-0.03999481201 


0.004740397135 


0.04620615641 


0.8677455584 


9 


0.01268447593 


0.009544584486 


-0.08312718427 


0.5234833748 


10 


-0.08279821256 


0.01060011122 


0.3628878114 


-5.029520772 


11 


0.02792425547 


0.02405914339 


-0.3008789073 


5.515860626 


12 


-0.1571703581 


0.02491494020 


0.7014666549 


-7.989323468 


13 


0.05760043398 


0.05506860056 


-0.6772442574 


11.90025208 


14 


-0.3147585986 


0.06230568066 


1.703826430 


-25.49575660 


15 


0.1235842297 


0.1303351125 


-1.691767724 


36.94182961 


16 


-0.6430290985 


0.1554514202 


3.999895232 


-70.54713230 


17 


0.2684262835 


0.3123609133 


-4.227999178 




18 


-1.343451682 


0.3887314574 


9.417754413 





The advantage of introducing this notion is that the number of different vertex 
structures is much smaller than the number of topologically inequivalent graphs. 
Except for factors contributed by the vertices themselves, we sum the (positive) 
weights of all graphs with the same vertex structure. This yields weights that do 
no more depend on the bare coupling constants. For instance, in this way, the 
coefficients of ( [104| ) become 

<i(A)= E C(/^i'/^2,...)-n(^ (A))'\ (109) 

where the b^^j^ do not depend on the coupling constants. 

Some numbers of vertex structures are listed in Tab. |^. For a given lattice, 
including dimension and topology, and given symmetry group 0(N), we don't need 
to keep all the graphs. Only those vertex structures together with the small set 
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Table 4: As in Tab. 3, but on the 6 x oo^ lattice. 



L 


iw 




iPl 


IW 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 


0.25 


-0.125 
0.01041666667 
0.01822916667 
0.0006510416667 
-0.03359646267 
0.007688395182 
-0.02836252848 
0.009167494597 
-0.05681311643 






0.0078125 


0.0004882812500 
0.002115885417 
0.005743408203 
0.004374186198 
0.007763120863 
0.008740446303 


-0.0625 


0.1875 
-0.04166666667 
-0.1562500000 
0.02148437500 
0.1898695204 
-0.07582329644 
0.004855290166 
-0.05144329777 
0.2625076545 


0.078125 


-0.6640625000 
0.2669270833 
1.263834635 
-0.3728027344 
-2.278199090 
1.235062069 
1.245049212 
0.06869500487 
-4.044370528 


11 
12 
13 

14 


0.01901101339 
-0.1080852959 
0.03790109259 
-0.2258664596 


0.01786831712 
0.01844869455 
0.03792749792 
0.04196657999 


-0.2090660888 
0.4873771492 
-0.4478781738 
1.260000420 


4.051910630 
-5.358685469 

7.732701464 
-18.95164927 


15 
16 
17 
18 


0.08196611346 
-0.4902317790 
0.1833679347 
-1.094178613 


0.08673188894 
0.09916992312 
0.2045461221 
0.2415480079 


-1.145127753 
3.140756286 

-2.979792099 
7.829428325 


25.74367205 
-55.14063511 



of numbers /i2, . . .) are needed for the final evaluation of the LCE series 

coefficients such as ( |109| ), for all values of the coupling constants A. 

As an example of the resulting series we consider the 0(4)-symmetric $^-theory, 
as described by the action ([l|)-@, with Xi = oo and A2 = 0. The coefficients of the 
various IPI susceptibilities are listed in the tables Tab. ^ and ^ for the hypercubic 
lattices of the sizes 4 x 00^, 6 x 00^ and 00'^, respectively. 

It appears that the IPI series are often of better use than the connected ones. 
For instance, the critical line ^^(Ai) can be determined as the smallest positive root 

°^ 1 - {2 K){2D)x r{Xi 

with significantly better accuracy than from the various ratio and Fade' methods 
applied to the series of X2(Ai). 

The four-dimensional finite temperature data yield the following numbers on the 
standard high temperature critical exponents: 



X2{x,r = ' (110) 
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Table 5: As in Tab. 3, but on the 3-dimensional oo^ lattice. 



L 


iPl 


f^i:2 


iPi 


iPl 





0.25 





-0.0625 


0.078125 


1 














2 


-0.09375 





0.140625 


-0.4980468750 


3 


0.0078125 


0.0078125 


-0.03125 


0.2001953125 


4 


0.01367187500 





-0.09472656250 


0.7391357422 


5 


-0.0004882812500 


-0.0004882812500 


0.01708984375 


-0.2545776367 


6 


-0.008585611979 


0.001464843750 


0.06603190104 


-0.8673782349 


7 


0.001572672526 


0.001572672526 


-0.02043457031 


0.3948656718 


8 


0.0003199259440 


0.001302083333 


-0.02937910292 


0.6583463351 


9 


0.0001472897000 


0.0005287594265 


0.007233217027 


-0.3005849621 


10 


-0.001710902320 


0.0008281707764 


0.02126218503 


-0.5705569410 


11 


0.0005031175714 


0.0009918579980 


-0.01068644069 


0.3722756644 


12 


-0.0007259439854 


0.0008832578306 


-0.002698613418 


0.2289861101 


13 


0.0002879108938 


0.0008595423960 


-0.001406942080 


-0.1410064280 


14 


-0.0007758281759 


0.0008468480700 


0.007465842599 


-0.1955813214 


15 


0.0002983264770 


0.0008996308209 


-0.004516607350 


0.2479154031 


16 


-0.0005740168913 


0.0008895557789 


0.005824725792 


-0.3629113447 


17 


0.0002457315931 


0.0008854585355 


-0.01304684165 




18 


-0.0006825215700 


0.0008996335601 


0.003926613713 





7 = 1.443(25), 

ly = 0.735(11), (111) 
r] = 0.035(28), 

whereas the three-dimensional series yield 

7 = 1.449(16), 

jy = 0.735(8), (112) 
77 = 0.029(17). 

Derivation and discussion of the above exponents as well as other critical data and 
the presentation of the methods used will be given elsewhere. Our aim here was to 
present those techniques that are capable of pushing the linked cluster expansion 
series towards very high orders, in such a way that details on critical behaviour can 
be obtained, even on asymmetric (finite temperature) lattices. 
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6 Extended vertex ordering for non-exact ring 
graphs. 

In this appendix we hst the notions needed to define the order relation O of ( ^7D for 
the vertices of a non-exact ring graph, as announced in subsection p.3| . Essentially, 
the notion of an internal 2- vertex is to be replaced by symmetric 2- vertex. 
Hence let P be a non-exact ring graph. We define 

O =< Olw, [Ons-, Osymm] > ■ (113) 
(2) 

Osymm is an order relation on B^ gy^^, the set of symmetric 2-vertices of P, andC„s 
an order relation on By \ -Bp^sj/mm- They are defined along the following lines. 

A symmetric 2-chain of f G -Bp^sj/mm is a path p = {wi, W2, ■ ■ ■ , Wn) with wi = v 

(2) 

and W2, . . . , Wn-i G B]. symm niutually distinct vertices. If the last vertex is not a 

(2) 

symmetric 2-vertex, Wn ^ B^ ^y^j^, p is called maximal symmetric 2-chain of v. If 

(2) 

V G Bj. gyj^j^, we assign a tripel of two integers and a vertex to p by 

M = n-l, (114) 

A(p) = c{v), and (115) 

e{p) = Wn- (116) 

(2) 

Every v G BY ^ymm has precisely two distinct maximal symmetric 2-chains. We 
denote them by Ci(f) and C2(f). They intersect in v and possibly the last vertex 
only. Furthermore we set 

(J^symmiv) = maX /i(Q). (117) 

Let be the order relation on Br \ B^^^^^^, defined by =< O5, Ons >, 

where 

v>w <^==^ max k > max k, (118) 

for all f , w G i3r \ Brlymm- 1^ induces an order on the maximal symmetric 2-chains 
{ci{v),C2{v)} by 

fi{ci) < fi{cj) , or 

Ci{v) > cj{v) /^(ci) = Kcj) and A(ci) > \{cj) , or (^g) 

/i(ci) = /i(cj) and A(ci) = A(cj) and e(ci) >"e(cj). 



34 



If unique, we choose "short and long chains" by 

Cs(v) = max cAv), 

1=1,2 

ci{v) = mill Ci{v). (120) 

If neither Ci{v) > C2{v) nor C2{v) > Ci{v), we choose any one of the maximal 
symmetric 2-chains as Cs{v) and the other one as ci{v). Again, this freedom is 
allowed because all that will be used is the order relation ( |119| ). 

C^o„ induces through ( p^ an order relation Osymm on B^^^^^^ by Osymm =< 
Oii,Oi2 >, with 

v>w <^=^ Cs{v) > Cs{w), (121) 
v>w ci{v)>ci{w), (122) 

for all symmetric 2-vertices v,w. Finally, Ons is defined as follows. With 

■5l(^) = J2 f^symmiu;), (123) 

i ,symm 

S2{v) = max Hsymm{w), (124) 

i ,symm 

(2) 

for every v ^ ^ , we set Ons =< C„s,i, Cns,2 >, where 



V y-'^'w Skiv) > Skiw), (125) 



for k = 1,2 and for all v,w &By \ B^^^y^^. 
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